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PREFACE 


The  U.S.  Air  Force  Project  RAND  research  program  at  The  Rand  Cor- 
poration has  for  many  years  included  studies  of  both  target  acquisition 
problems  and  weather  effects  on  military  operations.  The  work  described 
in  this  report  was  supported  by  Project  RAND  under  the  study  project 
"Science  and  Technology  Research."  It  is  a continuation  of  the  weather- 
effects  studies  with  application  to  target-acquisition  modeling. 

In  the  mainstream  of  target-acquisition  modeling — at  Rand  and  else- 
where— weather  factors  have  traditionally  entered  in  very  simple  terms 
that  are  difficult  to  relate  to  real,  dynamic  weather  situations  replete 
with  significant  time  and  space  variability.  To  help  correct  this  dif- 
ficulty, the  author  has  developed  a method  by  which  the  quality  of  the 
atmospheric  transmission  of  visible  and  infrared  radiation  (and  imagery) 
can  be  quantitatively  estimated  directly  from  the  hour-by-hour  weather 
observations  that  are  electronically  archived  for  thousands  of  locations 
worldwide.  Coupled  with  a target  acquisition  model,  sensor  characteris- 
tics, and  tactical  scenario,  the  method  peirmits  the  direct  evaluation 
of. weapon  systei  '•rformance  in  any  weather  situation  or  climatic  regime 
that  can  be  extracted  as  a subset  of  the  weather  data  base.  One  such 
evaluation  recently  completed  at  Rand  using  this  modeling  approach  ir 
reported  in  a forthcoming  Rand  report  on  the  utility  of  an  adverse- 
weather  precision-guided  munition  in  a NATO  context. 

As  indicated  by  its  subtitle  ("Weather  and  Warplanes  VI") , this 
report  is  preceded  by  five  others  that  deal  with  the  effects  of  weather 
and  weather  information  on  military  systems  and  operations: 


R-740-PR,  Use  of  Weather  Information  in  Determining  Cost/ 
Performance  and  Force-Mix  Tradeoffs:  Weather  and  War- 

planes If  R.  E.  Huschke,  June  1971  (Unclassified). 

R-742-PR,  Ten  Guidelines  for  the  Sirndation  of  Weather  Sen- 
sitive Military  Operations:  Weather  and  Warplanes  II, 

R.  E.  Huschke,  June  1971  (Unclassified). 


R-774-PR,  A Simple  Model  to  Elucidate  the  Utility  of  Weather 
Forecasting  in  Military  Operations:  Weather  and  Warplanes 

III,  R.  R.  Rapp,  August  1971  (Unclassified). 
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R-1195-1-PR,  Tactical  Ai^power  in  NATO  Contingencies — Modeling 
Weather  Constraints  on  Air  Operations:  Weather  and  War- 

planes IV  (U) , R.  ■'..  Huschke,  August  1973  (Confidential). 


R-1349-PR,  Main  Conclusions  from  the  'Weather-85'  Study: 
Weather  and  Warplanes  V (U) , R.  E.  Huschke,  January 
1974  (Secret). 


This  work  should  be  of  interest  to  Air  Force  and  other  DoD  agencies 
concerned  with  assessing  the  effects  of  weather  conditions  on  visual 
and  8-12  um  infrared  sensor  systems. 
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SUMMARY 


New  and  remarkably  effective  weapons  use  visible  and  infrared 
imaging  sensors  to  locate  targets  and  automatically  guide  the  missile^ 
to  their  targets  with  great  precision.  The  trouble  is  that  weather  con- 
ditions can,  and  in  some  parts  of  the  world  often  do,  deny  these  sys- 
tems their  needed  capability  to  see  through  the  atmosphere.  A problem 
that  arises,  then,  is  how  to  predict  the  performance  of  different  sen- 
sor systems  in  the  weather  conditions  that  characterize  different  places 
and  times — that  is,  in  the  climates  of  a variety  of  potential  theaters 
of  battle. 

World  weather  and  climate  are  depicted  in  most  relevant  temporal 
and  spatial  detail  by  huge  archives  of  surface  weather  observations 
from  which  subsets  can  be  taken  to  represent  almost  any  desired  region, 
season,  or  time  of  day.  From  the  information  contained  in  surface 
weather  observations,  it  is  possible  to  infer  visible  and  8-12  gm  image 
transmission  parameters,  thereby  opening  the  door  to  predicting  sensor 
performance  at  (almost)  any  place  of  interest,  including  seasonal  and 
diurnal  effects. 

Such  common  observables  as  cloud  amount  and  height,  visibility, 
precipitation,  wind  speed,  temperature,  dewpoint,  and  relative  humidity 
permit  construction  of  algorithms  to  make  quantitative  estimates  of  the 
following: 

o Visible  extinction  coefficient. 

o Sky-ground  ratio. 

o 8-12  gm  extinction  coefficient  (including  the  aerosol  extinc- 
tion coefficient  as  well  as  the  water  vapor  absorption  coeffi- 
cients) . 

o Atmospheric  vertical  profiles  of  the  relevant  variables. 

Some  of  these  algorithms,  notably  the  ones  for  sky-ground  ratio 
and  8-12  gm  aerosol  extinction,  warrant  strengthening  by  more  field 
measurement  studies  and  a broader  application  of  theory  than  has  been 
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INTRODUCTION 


A trend  in  the  U.S.  and  NATO  military  forces  is  to  rely  on  effec- 
tive weapons  to  help  counteract  the  numerical  superiority  in  armor  and 
manpower  held  by  Warsaw  Pact  forces.  A major  class  of  these  new,  ef- 
fective weapons  is  the  precision  guided  munitions  (PGMs)  that  use  optical 
sensing  devices — the  eye,  television,  or  imaging  Infrared  (HR) — for 
target  acquisition,  lock-on,  and  terminal  guidance.  Such  weapons  vary 
in  performance  because  weather  varies  in  its  degree  of  interference 
with  the  signals  as  they  propagate  from  target  to  sensor.  If  a target 
can  be  effectively  "seen"  through  the  acquire/lock-on/launch/guide  se- 
quence, expected  performance  is  excellent;  otherwise,  performance  is 
nil. 

Modern  tactical  warfare  is  apt  to  be  more  weather-sensitive  than 
in  the  past,  even  as  recently  as  in  Vietnam.  Weather  has  become  a fac- 
tor to  be  reckoned  with  seriously  at  all  levels  of  weaponry  decision- 
making. Procurement  and  deployment  decisions  require  statistically 
valid  predictions  of  PGM  utility  as  a function  of  weather  conditions 

f 

in  different  potential  theaters  of  combat  at  different  times  of  year. 
Day-to-day  and  hour-to-hour  d<  lent  decisions  require  adaptation  of 
immediate  weather  forecasts  ro  weapon  loadout  and  target  selection  de- 
cisions. 

The  time  and  space  variability  of  weather  is  well  documented  in 
a voluminous  data  base  containing  many  years’  history  of  hour-by-hour 
surface  weather  observations  taken  at  thousands  of  locations  worldwide. 
Each  observation  depicts  the  local  state  of  the  atmosphere  but  does  not 
explicitly  include  several  quantities  that  are  needed  to  julge  atmo- 
spheric effects  on  visual  and  infrared  target-acqulstion  and  missile- 
guidance  systems.  If  these  missing  quantities  could  be  inferred  from 
the  data  that  are  contained  in  the  observations,  then  the  surface 
weather  data  base  could  become  a valuable  resource  in  making  statisti- 
cally valid  predictions  of  weapon  performance  (including  diurnal,  sea- 
sonal, and  geographic  variability). 

The  model  of  weather  e^ffects  on  ^actical  J^arget  Requisition  (WETTA) 


tv 
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described  in  this  report  is  a collection  of  algorithms  by  which  it  is 
possible  to  estimate  quantities  used  in  mathematical  formulations  of 
visible  and  infrared  radiation  and  image  (contrast)  transmission  through 
the  atmosphere  from  surface  weather  data.  Some  of  the  algorithms  must 
be  regarded  as  first  approximations  because  they  are  attempts  to  quan- 
tify complex  phenomena  about  which  neither  empirical  nor  theoretical 
understanding  is  yet  sufficient  to  support  firmly  confident  quantifica- 
tion. 

Appendix  A details  the  portions  of  the  Balley-Mundie  [1,2]  visual 
and  infrared  target  acquisition  models  that  have  been  adapted  for  uj:e 
with  the  WETTA  algorithms  and  applied  to  several  practl'*al  probJems, 
including  a comparison  of  the  utilities  of  TV  (E-0)  and  iiraging  infrared 
(HR)  Maverick  missiles  in  the  context  of  an  air-ground  Ijattie  simula- 
tion used  on  a NATO  central  front  scenario  (3) . 
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II.  SURFACE  WEATHER  C3SERVATI0NS 


The  major  weather  seirvices  in  the  United  States  collect,  use,  and 
electronically  file  surface  weather  observations  (as  well  as  other 
types  of  observational  data)  originating  throughout  the  world.  For 
most  defense-related  applications,  the  primary  source  of  these  data  is 
the  United  States  Air  Force  Environmental  Technical  Applications  Center 
(USAFETAC) , Scott  Air  Force  Base,  Illinois,  operated  by  the  USAF  Air 
Weather  Service. 

Over  a period  of  about  seven  years.  The  Rand  Corporation  has  ac- 
cumulated a subset  of  the  surface  weather  data  base  from  USAFETAC  for 
specific  application  in  Rand  studies.  The  Rand  Weather  Data  Bank 
(RAWDAB)  files  now  contain  surface  weather  data  on  93  locations  in  22 
countries;  the  observation  Intervals  are  mostly  one  hour  and  three  hours, 
and  the  average  period  of  record  is  about  12  years  [4].  Table  1 lists. 


Table  1 


RAWpAB  FILES  FOR  GERMANY 


Location 

Dates  of  Record 
(mo/yr) 

Observation 

Interval 

(hr) 

Number  of 
Observations 

From 

To 

Berlin 

3/46 

12/63 

1 

155,468 

Bitburg 

3/52 

12/67 

1 

135,521 

Bremerhaven 

1/49 

11/71 

3 

45,424 

Dresden 

1/52 

12/63 

3 

14,298 

Emden-Hafen 

4/60 

11/71 

12,902 

Erfurt 

1/59 

12/63 

3 

12,535 

Essen 

1/49 

11/71 

3_ 

62,268 

Fulda 

9/60 

12/70 

1® 

58,854 

Grafenwohr 

6/62 

12/70 

1 

75,247 

Hamburg 

1/49 

11/71 

1,  3 

120,204 

Hannover 

1/49 

11/71 

1,  3 

94,094 

Heidelberg 

4/51 

12/70 

1 

172,595 

Kitzingen 

7/63 

12/70 

1 

62,252 

Leinefelde 

1/52 

12/60 

6 

13,025 

Leipzig 

1/52 

12/60 

6 

12,907 

Magdeburg 

1/52 

12/63 

3,  6 

22,796 

Munster 

8/59 

11/71 

3 

32,594 

Neublberg 

2/46 

1/58 

1 

104,778 

^ai.ily  daytime  observations. 


for  illustration,  the  RAWDAB  files  for  West  and  East  Germany,  including 
the  dates  of  record,  observation  interval,  and  total  number  of  observa- 
tions for  each  location.  The  fact  that  locations  having  the  same  dates 
of  record  and  observation  Interval  have  different  numbers  of  observa- 
tions (e.g.,  Hamburg  and  Hannover)  serves  as  a warning  that  these  data 
sets  commonly  contain  both  sporadic  and  systematic  gaps  that  the  user 
should  be  aware  of. 

Table  2 lists  the  meteorological  information  that  is  normally 
found  in  a surface  weather  observation.  Those  marked  with  "a"  are  the 
variables  on  which  the  visible  and  IR  transmission  algorithms  are  based. 
Celling  height  may  or  may  not  be  included  and  cloud  layer  data  may  be 
more  or  less  detailed,  both  depending  on  the  type  of  weather  code  for- 
mat used  in  the  original  compilation  and  transmission  of  the  observa- 
tion. In  the  present  work,  it  Is  assumed  that  ceiling  height  is  either 
given  or  is  deduced  from  available  cloud  data. 


Table  2 

INFORMATION  NORMALLY  CONTAINED  IN  A SURFACE 
WEATHER  OBSERVATION 

Atmospheric  State  Variables 

Atmospheric  pressure  at  station  elevation 
Atmospheric  pressure  reduced  to  sea  level 
Character  and  amount  of  3-hcur  pressure  change 
Temperature  (dry  bulb)® 

Temperature  (wet  bulb) 

Dewpoint  temperature® 

Relative  humidity^ 

Wind  direction 
Wind  speed® 

Weather  Phenomena  and  Obstructions  to'  Vision 
Visibility^ 

Present  weather  (includes  types  and  intensities  of 
obstructions  to  vision  and  precipitation)® 

Cloud  Information 

Total  cloud  cover® 

Ceiling  height  (height  of  that  cloud  layer  above 
which  less  than  half  the  sky  is  observable)® 

Cloud  layer  amounts® 

Cloud  layer  heights® 

Cloud  layer  types 

^ata  used  in  WETTA  algorithms. 
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III.  WETTA  MODEL  ALGORITHMS 


There  are  a number  of  target  acquisition  models  for  visual  and  in- 
frared systems  [1,2, 5-9],  and,  although  they  differ  somewhat  in  form 
and  complexity,  they  all  have  to  contend  (at  least  implicitly)  with  the 
same  sets  of  weather  conditions  and  phenomena  that  affect  the  transmis- 
sion of  image  information  from  target  scene  to  sensor.  I have  isolated 
these  universal  weather  factors  from  the  other  factors  involved  in  tar- 
get acquisition  modeling.  My  method  for  estimating  the  values  of  the 
relevant  atmospheric  parameters  is,  therefore,  independent  of  any  spe- 
cific target  acquisition  model. 

The  relationships  between  atmospheric  parameters  and  target  ac- 
quisition system  performance  are  complicated.  The  effects  of  the  dif- 
ferent atmospheric  parameters  are  interactive;  they  cannot  be  calculated 
separately  and  then  simply  combined,  either  for  a single  type  of  acqui- 
sition system  or  for  different  types  of  systems.  The  atmospheric  pa- 
rameters themselves  are  correlated  and,  hence,  cannot  be  treated  sta- 
tistically as  independent  variables.  Finally,  the  parameters  are 
highly  variable  in  space  and  time;  system  performance  is  apt  to  be 
similarly  variable,  which,  in  turn,  could  have  important  implications 
in  weapon-mix  deployments  and  in  weapon  use  concepts. 

The  main  problem,  then,  is  to  use  a source  of  weather  information 
in  which  as  much  as  possible  of  the  parametric  interdependency  and 
space-time  variability  is  intact.  The  computer-compatible  archives  of 
surface  weather  observations  meet  these  criteria  best.  In  essence,  the 
method  creates  a vertical-profile  "model"  of  atmospheric  parameters 
from  the  information  contained  in  each  surface  weather  observation. 

These  parameters  have  been  put  into  a target  acquisition  model  to  quan- 
tify the  performance  of  specific  systems  against  specific  targets  in  the 
climates  depicted  by  selected  sequences  of  historical  weather  data. 

In  developing  the  algorithms,  1 have  focused  on  the  visible  and 
8-12  pm  infrared  portions  of  the  electromagnetic  spectrum,  for  these 
are  the  principal  wavelength  regions  currently  being  exploited  for 
imaging  sensors.  The  image- transmission  related  parameters  estimated 
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from  surface  weather  observables  are  grouped  according  to  spectral  re- 
gion. 

Visible  Image  Transmission 

1.  Visible  extinction  coefficient 

2.  Sky-ground  ratio 

8-12  pm  Infrared  Image  Transmission 

3.  8-12  pm  extinction  coefficient 

a.  Water  vapor  molecular  absorption  coefficient 

b.  Water  vapor  continuum  absorption  coefficient 

c.  Aerosol  extinction  coefficient 
Both  Visible  and  Infrared  Transmission 

4.  Cloud-free  line  of  sight 

5.  Atmospheric  vertical  structure 


VISIBLE  EXTINCTION  COEFFICIENT 

"Visibility"  is  the  routinely  observed  surrogate  for  the  visible 
extinction  coefficient.  Visible  transmission  is  affected  by  all  atmo- 
spheric aerosols  (haze,  smoke,  smog,  fog,  etc.).  The  Koschmieder  visi- 
bility theory  (10]  defines  the  relationship, 

where  V is  the  maximum  range  at  which  an  ideal  black  target  can  be  seen 
against  the  horizon  (the  "visibility"),  B is  the  visible  extinction  co- 
efficient, and  C^  is  the  observer's  threshold  contrast  (the  contrast 
ac  which  50  percent  of  all  observers  would  see  the  target) . Middleton 
[11]  discusses  at  length  the  practical  problems  in  using  this  equation 
to  estimate  3 from  observations  of  V.  Visibility  is  estimated  by  ob- 
servers having  different  threshold  contrasts  using  nonstandard  arrays 
of  ncr.i-deal  visibility  targets. 

Koschmieder  assumed  = 0.02,  whence 


V = 


3.912 

3 
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whlch  has  become  adopted  as  a sort  of  "standard"  and  been  given  various 
names,  such  as  "meterological  range,"  to  distinguish  it  from  observed 
visibility.  Measurements  by  Hering  [12]  stiongly  support  the  ■ 0.02 
value  for  ideal  visibility  estimating  conditions  (a  carefully  placed 
array  of  large  black  markers) , but  the  effective  value  of  increased 
rapidly  when  natural  markers  beyond  the  array  were  used.  Middleton  re- 
ports two  comparisons  of  contrast  measurements  against  visibility  ob- 
servations under  field  conditions,  resulting  in  mean  effective  threshold 
contrasts  of  0.031  and  0.039.  This,  plus  a smattering  of  additional 
evidence  also  reported  by  Middleton,  leads  me  to  adopt  the  value  = 
0.035  as  a statistically  representative  bridge  between  reported  visi- 
bilities and  actual  extinction  coefficients;  therefore,  for  present 
purposes , 


8 


3.352 

V 


If  V is  bm,  the  units  for  8 are  neper  km  ^ (the  neper  is  the  natural 
logarithmic  analog  of  the  decibel).  The  difference  between  these  two 
assumptions  for  is  shown  in  Fig.  1.  The  use  of  = 0.035  results 
in  the  inferred  visible  extinction  coefficient  being  15  percent  lower 
than  with  = 0.02. 

SKY-GROUND  RATIO 

Visual  image  transmission  is  formulated  as  target-to-background 

contrast  transmission,  per  Duntley  [13].  His  algebraic  formulation 

makes  use  of  the  ratio  of  the  luminance  of  the  horizon  sky  (L^)  to  the 

luminance  of  the  ground  near  the  target  (L  ) , both  taken  in  the  azimuth 

8 1 

toward  the  target,  the  so-called  "sky-ground  ratio."  It  arises,  phys- 
ically, because  of  extraneous  light  being  scattered  into  the  sensor 
field  of  view  by  atmospheric  aerosols.  The  received  contrast  at  range 


Some  researchers  prefer  the  use  of  a closely  related  variable, 
"directional  path  reflectance,"  in  lieu  of  the  "sky-ground  ratio." 

The  preference,  in  my  opinion,  is  arbitrary,  for  both  involve  similar 
problems  in  measurement  or  approximation.  For  backgrounds  and  targets 
whose  luminances  consist  entirely  of  reflected  "'ight,  their  reflectances 
can  be  used  in  place  of  their  luminances;  and  we  do  so  commonly  when 
dealing  with  militarv  target  detection  in  daylight. 
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sun,  observer,  and  target),  extinction  coefficient  (visibility),  cloud 
conditions,  and  surface  albedo  (which  I take  to  be  the  same  as  back- 
ground luminance).  There  is  no  rigorous  analytical  way  to  calculate 
Lg/L^  from  weather  data;  therefore,  a purely  empirical  method  is  used, 
based  on  calculations  given  in  Duff  [14]  who  used  both  radiometric  mea- 
surements (by  Duntley  et  al.  [15]  in  Southern  Germany)  and  the  RRA/ 

AFCRL  model  of  multiple  scattering  in  the  atmosphere  [16].  A detailed 
examination  of  Duff's  calculations  permits  the  following  generalizations 
to  be  drawi : The  sky-ground  ratio 


(a)  increases  with  decreasing  surface  albedo; 

(b)  increases  with  decreasing  visibility; 

(c)  increases  with  decreasing  solar  elevation  ajigle; 

(d)  is  independent  of  sun  angle  under  cloudy  skies; 

(e)  is  maximum  with  receiver  depression  angle  (look  angle) 

10-30“; 

(f)  is  minimum  witli  receiver  depression  .ingle  near  0°  and 
50-70“; 

(g)  has  maximum  value  about  2 to  3 times  minimum  values  as 
a function  of  receiver  depression  angle; 

(h)  does  not  aepend  strongly  on  azimuth  relative  to  the  sun; 
and 

(i)  averages  about  2 to  5 in  Germany  in  summer  with  good 
visibility. 

The  sky-ground  ratio  algorithm  produces  results  in  conformity  with 
the  above  general  rules.  First,  an  estimate  for  sky-ground  ratio  (SG) 
is  made  independently  for  each  of  three  of  the  four  variables  of  which 
it  is  a function.  These  three  are  the  surface  albedo  (p) , the  solar 
elevation  angle  (c)  and  the  visibility  (V).  In  each  case,  the  shape 
or  trend  of  the  functional  relationship  as  calculated  by  Duff  is  not 
violated,  and  values  are  reproduced  that  are  consistent  with  those  that 
correspond  to  a receiver  depression  angle  (6),  the  fourth  variable,  be- 
tween about  80“  and  50“ . These  independent  estimates  are  averaged  and 
then  multiplied  by  a factor  (F^)  that  accounts  for  the  dependence  of 
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L /L  on  6,  the  form  of  which  varies  as  a function  of  the  first  three 
variables.  Hence,  a sky-ground  ratio  estimate,  SG,  is  calculated  as 
follows ; 

L /L  « ^ - F-  + ^u)/3, 

s g 0 p e V 

where 

^ = exp  (-1.15  M p - .75), 

P 

SG  = exp  (-3.4  sin  e + 2.7), 

= exp  (-.5  V 4-  L.5). 

If  the  sky  is  broken  or  overcast  (cloud  amount  > 4/8),  e is  set  at  25°. 
Curves  for  the  three  expressions  for  SG  are  plotted  as  solid  curves  in 
Fig.  2.  The  dashed  curves  on  Figs.  2a  and  2b  are  Duff's  calculations 
from  the  RRA/AFCRL  multiple-scattering  model  for  values  of  the  otlier 
variables  that  more  or  less  bound  my  approximation.  The  curve  with 
open  circles  on  Fig.  2a  is  a long-standing  approximation  (due  to  Duntley 
as  reported  by  Middleton  [11,])  of  sky-ground  ratio  as  a function  of 
surface  albedo;  ^ is  a compromise  between  Duff's  calculations  and 
Duntley 's  suggested  approximation.  As  far  as  the  visibility  function 


Fig.  2 — Independent  estimatej  of  jky-ground  ratio  (SG)  from 
(a)svirfare  albedo,  (b)  solar  elevotion  angle,  and  (c)  visibility 
with  other  estimates  and  colculations  shown  for  reference  purposes. 
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Is  concerned  (Fig.  2c),  very  little  infoinnation  was  available  on  which 
to  base  any  functional  relationship,  because  Duff's  calculations  were 
for  only  two  visibilities,  3 km  and  24  km,  resulting  in  ranges  of  values 
fo.  sky-ground  ratio  as  shown  by  the  vertical  bars. 

The  factor  is  calculated  as  in  Tables  3 and  4,  which  were  con- 
structed empirically  from  examination  of  many  plots  (from  Duff's  calcu- 
lations) of  L /L  vs  6 for  different  combinations  of  e,  p and  V.  To 
9 8 

use  these  tables,  first  determine  values  of  indices  j^,  ^2*  ^3* 
from  the  assumed  or  calculated  values  of  6,  e,  p,  and  V as  in  Table  3. 
Then,  in  Table  4,  the  appropriate  value  of  F^  is  fcund  opposite  the 
resulting  combination  of  j-index  values.  Figure  3 shows  sample  plots 
of  sky-ground  ratio  (as  calculated  by  the  multiple-scattering  model) 
versus  receiver  depression  angle,  6,  for .two  visibilities,  two  albedos, 
two  sun  angles,  and  overcast  skies.  Also  plotted  are  data  points,  dots 
(averages) , and  vertical  bars  (ranges  of  values)  for  four  measurement 
flights  over  Southern  Germany,  for  which  p w .1,  35“  < e < 50“ , and 
7 km  < V < 30  km.  The  horizontai  dotted  lines  on  Fig.  3 represent  al- 
gorlthm  calculations,  SG,  that  correspond  to  values  of  the  variables  for 
which  the  uppermost  dashed  curve  was  calculated  by  the  multiple-scatter- 
ing model  (p  = .1,  e = 20*,  V = 25  km).  The  difference  between  the 
two  estimates  is  mainly  because  SG^  gives,  at  small  values  of  p,  a 
sky-ground  ratio  estlm£.ce  that  is  greater  than  that  estimated  by  the 
multiple-scattering  model.  Finally,  the  open  circles  on  the  vertical 
bars  are  algorithm  calculations  for  the  approximate  average  of  condi- 
tions encountered  by  the  measurement  flights  over  Southern  Germany 
(p  = .1,  e = 43“,  V = 18  km). 


8-12  urn  EXTINCTION  COEFFICIENT 

The  total  atmospheric  extinction  coefficient  for  8-12  pm  radiation 
can  be  taken,  to  a good  approximation,  to  be  the  sum  of  a water  vapor 
molecular  absorption  coefficient  (y^^^) , a water  vapor  continuum  absorp- 
tion coefficient  (y  ),  and  an  aerosol  extinction  coefficient  (y  ).  The 

C 3l 

first  two  are  functions  of  absolute  humidity;  the  second  is  also  a func- 
tion of  temperature. 
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At  greater  heights  up  to  «6000  m,  the  absolute  humidity, 

A Rj  A'  (1  + .0003  Z) , 

where  A'  is  calculated  from  the  dewpoint  at  height  Z (m)‘  as  ^i.  it  were 
the  dewpoint  at  Z = 0 («1000  mb).^ 


Water  Vapor  Molecular  Absorportion  Coefficient  (8-12  urn) 

McClatchey  et  al.  [17]  give  graphical  methods  by  which  the  two 

types  of  water  vapor  (HO)  absorption  coefficients,  v and  v , can  be 

Z -me 

estimated  from  the  absolute  humidity.  The  basis  of  these  methods  has 
also  been  incorporated  by  Selby  and  McClatchey  into  a computer  program 
for  atmospheric  transmittance  over  the  wavelength  range  0.25  - 28.5  ym, 
the  latest  published  version  of  which  is  LOWTRAN  3 (18). 

A curve  fit  to  the  average  H^O  molecular  absorption  coefficient 
over  8-12  ym  calculated  by  the  method  of  Ref.  17  gives 

^°«io  ^°^io  ^ 

Water  Vapor  Continuum  Absorption  Coefficient  (8-12  ym) 

When  the  method  of  Ref.  17  is  used  to  calculate  H^O  continuum 
absorption  coefficients  for  8-12  ym  radiation,  over  consistent  combina- 
tions of  temperatures  from  22  to  94®F  and  pressures  from  1000  to  800  mb, 
a reasonable  linear  fit  is  given  by 


'c(Ref.  17) 


.01432A. 


Roberts,  Blberman,  and  Selby  [19]  have  suggested  corrections  to 
the  Ref.  17  method  for  based  upon  their  examination  of  a large 
number  of  more  recent  atmospheric  and  laboratory  measurements.  First, 
they  suggest  an  approximate  25  percent  reduction  in  the  8-12  ym  cal- 
culated by  the  older  method  and  state  that  the  corrected  value  should 


Approximations  for  absolute  humidity  derived  from  basic  hygrome- 
trlc  formulas.  A'  is  accurate  within  5 percent  in  the  dewpoint  range 
10°F  - 70“F. 
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be  valid  at  296°K  (73®F).  They  further  suggest  a strong  temperature 
dependence  of  y^.  (These  corrections  have  been  accepted  by  the  de- 
signers of  the  computer  model  LOWTRAN  3,  and  are  incorporated  into  a 
revision,  LOWTRAN  3A) . The  resulting  approximate  value  is  given  by 


« .01074A  exp  (1800 


T + 255.4 


1 

296 


where  T is  temperature  in  °F. 

A graphical  summary  of  these  approximations  for  both  types  of 

H2O  absorption  coefficients  is  presented  in  Fig.  4.  "X"  plots  and  the 

vertical  bars  indicate,  respectively,  values  calculated  for  y and 

m 

ranj.es  of  values  calculated  for  y^  by  the  graphical  method  of  Ref.  17; 

and  the  associated  curves  are  produced  by  approximating  expressions. 

The  solid  curves  for  y incorporate  the  corrections  of  Ref.  19,  and  are 

c 

terminated  at  the  point  where  the  combinations  of  absolute  humidity  and 
temperature  are  equivalent  to  100  percent  relative  humidity.  The  cor- 
rections of  Ref.  19  are  incorporated  exactly  as  suggested  by  the  authors 
of  that  report;  therefore,  any  discrepancy  between  my  approximations 
and  LOWTRAN  3A  calculations  should  be  due  to  my  original  approximations 
to  the  Ref.  17  method. 

Aerosol  Kxtinction  Coefficient  (8-12  urn) 

Infrared  transmission  is  attenuated  by  atmospheric  aerosols  through 
botli  scattering  and  absorption.  Aerosol  extinction  (scattering  plus 
absorption)  is  a strong  function  of  particle  size  relative  to  the  wave- 
length of  the  radiation,  with  significant  extinction  occurring,  in  gen- 
eral, only  if  particle  radii  are  of  the  same  order  of  magnitude  as,  or 
larger  than,  the  wavelength.  As  a result,  the  IR  aerosol  e;  .inccion 
coeffirieni  is  about  two  orders  of  magnitude  smaller  than  the  visible 
extinction  coefficient  in  small-particle  aerosols  (such  as  a dry  con- 
tinental haze),  but  is  about  equal  to  the  visible  extinction  coefficient 
in  large-particle  aerosols  (such  as  heavy  fog,  clouds,  and  precipita- 
tion) . 


/c  / 


(Ref.  17) /■ 

X 

/■  ^r,(73^ 


7c  (53°F) 


X (23°F) 


(13°F), 
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Absolute  humidity  (gm  m~3) 

“ig.  4 — Approximations  for  H2O  molecular  (y^^)  and  continuum  (/j.)  cbsorption 

coefficients  for  8-12  radiation 


A very  simple  algorithm  has  been  developed  to  approximate  the  com- 
plex effects  of  atmospheric  aerosols  on  8-12  pm  transmission.  Aerosol 
particle  size,  a major  governing  factor,  is  a strong  function  of  rela- 
tive humidity  (RH) ; and  the  rate  at  which  a particle  grows  with  increas- 
ing Rll  depends,  further,  on  its  chemical  composition  and  initial  size 
(see,  e.g.,  Orr  ct  al.  [20];  Koenig  [21];  Hanel  [22]).  In  general, 
size  increases  per  incremental  increase  in  RH  are  small  with  RH  less 
than  about  60-75  percent,  and  the  crowth  rate  increases  in  an  exponen- 
tial manner  at  higher  Ril  to  about  louble  in  raeius  at  90-95  percent  RH 
and  about  quadruple  in  radius  at  RH  « 99  percent.  As  particles  grow 


w * 


with  increasing  RH,  the  aerosol  extinction  coefficient  at  8-12  wm_(7  ) 

di 

should  increase  more  rapidly  than  the  visible  extinction  coefficient 
(8);  also,  the  ratio  y /8  should  be  expressible,  to  a first  approxima- 
tion,  as  a function  of  RH. 

Several  investigators  have  shown  that  y^/3  « 1 in  fogs,  clouds  and 
precipitation,  all  of  which  are  phenomena  associated  with  RH  « 100  per- 
cent (see,  e.g.,  Deirmendjian  [23,24];  Ruppersberg,  bchellhase,  and 

Schuster  [251;  Rensch  and  Long  [26]).  Estimates  of  the  ratio  y /3  in 

a 

"dry"  atmospheric  aerosols  (RH  < 50  percent)  are  much  more  difficult 

to  find.  I have  assumed  that  Deirmendjian' s extinction  coefficient 
calculations  for  his  continental  "Haze  L"  particle  size  distribution 
model  are  representative,  in  general,  of  dry  atmospheric  aerosols. 

Based  thereon,  my  assumption  is  that  7 /8  w .015  for  RH  £ 50  percent. 

I have  further  assumed  that  Deirmendjian' s maritime  "Haze  M"  distribu- 
tion is  generally  representative  of  aerosols  with  RH  w 90  percent;  his 
"Haze  M"  calculations  give  y /3  ^ -065  [23].  A function  that  adequately 
fits  these  assumed  pairs  of  RH  and  y / values  is. 

Si 

for  50  percent  ^ RH  _<  99.5  percent, 


For  RH  < 50  percent. 


^ = .015; 


and  for  RH  > 99.5  percent,  or  when  precipitation  is  falling. 


~ E 1.0. 


This  set  of  functions  is  illustrated  in  Fig.  5. 


CLOUD-FREE  LINE  OF  SIGHT 

Table  5,  adapted  from  Rapp  and  Schutz  [27],  is  used  to  estimate 

the  probability  of  a cloud-free  line  of  sight,  P , as  a function 

Ci*  LOS 


P 0,1 


20  40  60  80 

Relative  humidity  (%) 

Fig.  5 — Assumed  function  relating  to  relative  humidity 


of  sensor  depression  angle,  6,  and  total  cloud  amount,  N(8ths),  in  all 

cloud  layers  between  tne  sensor  altitude  and  the  ground.  The  original 

data  Lor  are  those  of  Shanklin  and  Landwehr  [28],  which  Rapp 

Cr  LUb 

and  Schutz  aggregated  by  averaging  over  all  line-of-sight  azimuth  angles 
and  interpoiated  for  cloud  amounts  reported  in  8ths  rather  than  in  lOths, 
All  clouds  are  assumed  to  be  opaque  for  both  visible  and  8-12  ym 
radiation;  therefore,  for  practical  application  in  air-to-ground  target 
acquisition  problems,  their  pvinoipal  effect  is  to  limit  the  maximum 
altitude  of  air-to-ground  systems. 


ATMOSPHERIC  VERTICAL  STRUCTURE 

Aircraft  measurements  of  visible  extincticu  coefficients  in  the 
atmosphere  by  Duntley  et  al.  [11.29-31]  give  strong  support  to  the  view 
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Table  5 

PROBABILITY  OF  A CLOUD  FREE  LINE  OF  SIGHT  (PpFLOs) 
AS  A FUNCTION  OF  SENSOR  DEPRESSION  ANGLE  (6)  AND 
INTERVENING  CLOUD  AMOUNT  (N) 


N (8ths)  as 

Reported 

by 

Ground- 

-based  Weather 

Observers 

6^ 

0 

1 

2 

3 

4 

5 

6 

7 

8 

0^ 

(.94) 

(.78) 

(.63)  ( 

.48 

(.36) 

(.27) 

(.18) 

(.09) 

(0) 

10 

.96 

.85 

.71 

.59 

.46 

.36 

.27 

.15 

.01 

20 

.91 

.90 

.80 

.68 

.58 

.48 

.36 

.21 

.03 

30 

.98 

.92 

.85 

.76 

.66 

.55 

.43 

.26 

.05 

40 

.99 

.93 

.88 

.80 

.68 

.59 

.47 

.30 

.07 

50 

.99 

.94 

.89 

.82 

.74 

.64 

.50 

.33 

.08 

60 

1.00 

.95 

.90 

.83 

.76 

.67 

.53 

.34 

.08 

70 

1.00 

.95 

.90 

.84 

.77 

.68 

.54 

.35 

.08 

80 

1.00 

.95 

.90 

.84 

.77 

.69 

.54 

.35 

.08 

90 

1.00 

.95 

.90 

.84 

.77 

.69 

.55 

.35 

.08 

Linear  Interpolation  in  6 and  N is  permissible.  N will  nor- 
mally be  an  Integer  from  0 to  8. 

^Use  5=0  row  only  for  interpolating  between  5=0  and  5 = 
10.  For  a svirf ace-to-surface  look-path  (5  = 0),  P^p^Qj^  = 1.00. 


rhat  non-cloud  aerosols  (dust  and  other  particulate  pollutants)  are 
usually  homogeneously  distributed  throughout  a well-mixed  layer  of  the 
atmosphere  that  extends  from  the  earth's  surface  to  about  1500  m,  + 

1000  m,  approximately,  depending  on  the  gross  stability  characteristics 
of  the  atmosphere  at  the  given  time  and  place.  Of  26  extinction  coef- 
ficient profiles  measured  over  Southern  Germany,  Illinois,  and  Nev» 
.Mexico,  22  clearly  show  a homogeneous  boundary  layer  with  extinction 
coefficient  very  near  the  surface  value,  above  which  an  apparently  typi- 
cal upper-atmosphere  extinction  coefficient  of  0.02  - 0.06  neper  km'^ 
usually  prevails.  Because  of  this  and  other  considerations,  I have 
adopted  a layered,  rather  than  exponential,  model  of  vertical  atmo- 
spheric structure. 

Tlie  number  of  layers,  layer  dimensions,  and  physical  characteris- 
tics are  deduced  from  surface  weather  observables  by  the  decision  tree 
shown  in  Fig.  6.  The  determining  parameters,  all  strongly  associated 
with  the  static  stability  of  the  lower  atmosphere,  are  precipitation 


(whether  or  not  occurring),  cloud  amount,  wind  speed,  visibility,  and 
solar  elevation  angle.  Four  basic  vertical  structure  "models"  are 
used  in  which  the  atmosphere  is  divided  into  two  or  three  layers  of 
variable  thicknesses.  The  layer  properties  pertain  only  to  cloud-free 
portions  of  the  layers.  The  layers  are  horizontally  homogeneous,  and 
every  layer  is  vertically  homogeneous  in  its  aerosol  extinction  proper- 
ties. Temperature  and  dewpoint  lapse  rates  (rates  of  decrease  with 

altitude),  AT  and  AT,,  from  which  absolute  and  relative  humidities  are 

d 

calculated,  vary  according  to  atmospheric  stability  inferences.  The 
lowest  layer  always  has  the  surface  visibility  (visible  extinction  co- 
e^’icient),  and  the  temperature  and  dewpoint  profiles  are  always  "an- 
chored" by  their  surface  values.  The  following  are  brief  descriptions 
of  the  four  basic  structures: 

(a)  When  precipitation  is  falling,  the  lower  layer  extends  to 
the  base  of  the  precipitating  clouds.  An  upper  layer  is  defined  that 
pertains  to  the  cloud-free  air  between  and  above  the  precipitating 
clouds  (this  is  nearly  irrelevant  for  visible  and  IR  systems  since,  by 
fiat,  precipitation  from  scattered  clouds  is  not  treated  as  a "precipi- 
tation" case). 

(b)  Conditions  of  light  wind,  few  clouds,  and  low  (or  no)  sun  in- 
dicate a stable  surface  layer,  which  is  given  an  arbitrary  depth  of 
200  m.  The  single  upper  layer  has  excellent  visibility  and  is  made 
much  drier  than  the  lower  layer  by  imposing  a 10°C  dewpoint  disconti- 
nuity at  the  interlayer  boundary. 

(c)  When  above  conditions  are  not  met,  a neutral  or  unstable  at- 
mosphere is  indicated;  and  if  visibility  does  not  indicate  a heavy  fog 
or  local  heavy  pollution,  a lower  mixed  layer  is  assumed  with  depth 
equal  to  the  height  of  the  lowest  cloud  layer  (_<  2500  m and  > 200  m) 
or,  by  default,  1500  m.  The  upper  layer  has  excellent  visibility. 

There  is  nc  dewpoint  discontinuity  at  the  interlayer  boundary  as  in 
(b) , above. 

(d)  With  all  conditions  as  in  (c)  , above,  except  visibility  £0.5 
km,  it  is  assumed  that  there  is  a lower  layer  of  fog  or  pollution  (200 
m depth)  and  a middle  layer  with  fairly  good  visibility  (17  km).  The 
depth  of  the  middle  layer  and  all  characteristics  of  the  uppermost  layer 
are  calculated  as  in  the  two  layers  in  the  type  (c)  structure. 
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The  atmospheric  vertical  structure  models  thus  defined  for  each 
surface  weather  observation  contain  all  the  physical  variables  and 
spatial  dimensions  needed  to  calculate  the  atmospheric  transfer  func- 
tions of  visible  contrast  and  IR  signal-to-noise  ratio  along  any  spe- 
cified or  calculated  line  of  sight.  It  is  important  to  be  aware  that 
the  assumption  of  horizontal  homogeneity  and  the  simplistic  depiction 
of  elevated  layers  restricts  the  applicability  of  the  model  to  target 
acquisition  at  short  ranges  and  low  altitudes  (say,  < 20  km  range  and 
< 5000  m altitude) . 


I 
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IV.  NEEDS  FOR  IMPROVEMENT 


This  work  has  made  clear  the  need  for  improved  understanding  in  a 
number  of  areas: 

o Sky-ground  ratio y L /L  . Duff's  calculations  of  sky-ground 

s g 

ratio  [14]  are  unique.  Atmospheric  transfer  of  visual  con- 
trast appears  to  be  so  sensitive  to  this  parameter  (see  Fig. 
7)  that  more  measurements  and  theoretical  calculations  are 
warranted. 
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Fig.  7 — Calculated  visual  target  detection  probability  (Povis) 
as  a function  of  visibility  for  different  value  of  sky-ground 
ratio  (Ls/Lg).  Characteristic  target  dimension  is  3 m and 
inherent  contrast  is  0.4. 


o 8-12  vm  aerosol  '-xtinotion  aoeffioient  ^ y • Hy  scheme  to 
scale  this  parameter  as  a function  of  visible  extinction  co- 
efficient (visibility)  and  relative  humidity  is  tentative  and 
simplistic,  but  is  in  the  right  direction.  In  recent  tests 
of  HR  seekers  in  Germany,  serious  image  degradation  vas  con- 
sistently noted  coincident  with  low  visibility  and  high  rela- 
tive humidity  [32,33].  My  algorithm  says,  in  effect,  that 
dry  hazes,  smokes,  etc.  contain  mainly  small  particles  that 
degrade  visible  transmissions  much  more  than  IR  transmission, 
and  that  very  wet  aerosols  (RH  ->■  100  percent)  contain  enough 
large  particles  that  have  grown  through  accretion  of  water 
that  IR  degradation  approaches  visible  degradation.  I have 
set  aside,  at  least  for  the  present,  some  potentially  impor- 
tant differences  among  aerosols  of  different  origins  (sea  salt, 
continental  dust,  forest  fire  smoke,  industrial  effluents, 
etc.),  such  as  different  "dry"  size  distributions,  different 
growth  rates  with  relative  humidity,  and  different  spectral  ab- 
sorptivities.  In  the  context  of  my  modeling  approach — statis- 
tically valid  modeling  from  surface  weather  data — I do  not 
know  if  these  differences  would  be  significant.  Perhaps  each 
major  generic  type  of  aerosol  should  be  represented  by  its 
own  scaling  function  for  y ; the  decision  as  to  which  aerosol 
type  exists  would  then  have  to  be  based  on  geographical  and 
meteorological  criteria. 

There  are  also  miscellaneous  uncertainties  in  deducing  atmospheric 
vertical  structure  and  CFLOS  probability,  but  1 consider  these  to  be 
unimportant  compared  with  the  above  problems  and  compared  with  the  fol- 
lowing questions,  which  have  hardly  been  addressed  at  all  in  the  context 
of  target  detection: 

o Effects  of  aerosol  infrared  emissions.  Car Ion  [34,35]  has 
suggested  that  emissions  of  atmospheric  aerosols  at  infrared 
wavelengths  can  be  a major  factor  in  reducing  the  received 
contrast  of  thermal  images.  Absorption  by  liquid  water  is 
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strong  in  the  8-12  titn  region;  thermal  emission  by  the  absorber 
will  be  correspondingly  strong,  but  most  important,  the  wave- 
length of  maximum  emission  intensity  for  a black  body  at  typi- 
cal ambient  air  (and  aerosol)  temperatures  (r«  285°K)  is  near 
10  urn.  One  would  expect,  therefore,  that  when  a water  aerosol 
is  present  (nominally  at  BH  > « 60  percent,  especially  at 
RH  ->■  100  percent)  tne  atmosphere  would  "glow"  in  the  8-12  um 
band  and  have  a degrading  effect  on  IR  image  contrast  trans- 
mission that  is  similar  to  the  effect  of  scattered  "airlight" 
or  "path  luminance"  (parameterized  by  the  sky-ground  ratio) 
on  the  transmission  of  visual  image  contrast, 
o Wathev  effects  on  the  target  scene.  This  problem  encompasses 
both  visible  and  IR  depictions  of  the  target  scene,  but  I pass 
over  the  visible  problems  (shadow  variations,  snow-covered 
ground,  etc.)  in  favor  ol  the  problems  associated  with  deduc- 
ing the  thermal  characteristics  of  the  target  scene  from  sur- 
face weather  data.  Batten  [36]  has  suggested  a preliminary 
quantitative  algorithm  for  typical  temperature  differences 
found  among  the  natural  features  of  pastoral  terrain  in  north- 
temperate  latitudes,  based  on  season,  time  of  day,  and  weather 
conditions.  His  data  sources  were  two  sets  of  airborne  IR 
radiometer  measurements  (in  Wisconsin  [37]  and  in  Germany  [38]) 
and  one  set  of  measurements  made  with  tower-mounted  radiometers 
(in  New  Hampshire)  [39].  Batten's  algorithm  for  background 

thermal  clutter,  T,,  (°C),  is  summarized  in  Table  6 and  Fig. 

ab 

8.  None  of  these  basic  data  sets  were  acquired  with  HR  tar- 
get detection  problems  in  mind,  so  there  could  be  some  impor- 
tant interpretation  problems  relating  to  scene  content  and 
scale. 
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Table  6 

VALUES  OF  BACKGROUND  CLUTTER  FOR  VARIOUS  METEOROLOGICAL  CONDITIONS 


Meteorological  Condition 

"Ab 

Clear  to  scattered  skies 

Seasonal  variation  given  in  Fig,  8a  for 

(0  to  4/8) 

dawn  and  midday.  Diurnal  variation  for 

summer  and  winter  given  by  Fig.  8b. 

Broken  to  overcast  skies 

Use  dawn  values  given  in  Fig.  8a. 

(>  4/8  to  8/8) 

During  precipitation 

1°C 

The  24  hours  following 

3°C 

precipitation 

Snow  on  ground 

2“C 

Winds  greater  than  10  knots 

Use  dawn  values  given  in  Fig.  8a. 
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Fig.  8 Approximate  seosonal  ond  diurnal  voriations  of  bocVground 
Ibermol  clutter,  w for  pastoral  terroin  in  north*temperate 
lotitudei  with  cleor  sViet  and  low  winds 


V.  CONCLUDING  REMARKS 


The  foregoing  algorithms  were  developed  for  the  primary  purpose 
of  statistically  valid  prediction  of  visible  atid  8-12  urr.  transmissiot 
properties  of  the  atmosphere  in  a variety  of  climates  (different  loca- 
tions and  seasonal  and  diurnal  variations).  Computationally,  the  al- 
gorithms are  fairly  simple  because,  to  attain  the  goal  of  statistical 
validity,  it  is  necessary  to  perform  calculations  on  very  large  numbers 
of  past  surface  weather  observations.  To  date,  however,  the  model  has 
been  applied  only  to  selected  small  samples  of  data. 

The  methodology  described  in  this  report  clears  the  way  for  the 
application  of  historical  weather-data  sets  to  the  evaluation  of  multi- 
spectral  "seeing"  conditions  for  almost  any  location  in  the  world.  The 
algorithms  undoubtedly  could  be  improved.  It  is  gratifying  to  note, 
in  this  regard,  tliat  most  new  programs  for  measuring  atmospheric  optical 
and  image-transmission  properties  (e.g. , Moulton  [33],  Fenn  [40],  and 
Hubbard  [41])  include  the  collection  of  standard  surface  weather  observa- 
tions as  part  of  their  basic  data  sets  for  the  specific  purpose  of  cor- 
relating the  weather  observations  with  the  optical  measurements. 

The  methodology  could  also  be  expanded  to  include  other  portions 
of  the  electromagnetic  spectrum  in  which  sensors  do  or  might  operate. 
Inclusion  of  the  3-5  urn  Infrared  band  would  be  quite  simple.  Expansion 
of  the  methodology  to  the  near -millimeter  and  microwave  spectral  re- 
gions, however,  would  be  much  more  difficult.  At  those  longer  wave- 
lengths, clouds  cannot  be  assumed  to  be  opaque,  and  transmission  is  a 
strong  function  of  cloud  and  precipitation  particle  size  distributions 
and  concentrations.  Modeling  of  these  parameters  from  surface  weather 
observations  might  be  worth  attempting,  but  the  uncertainties  would  be 
large. 
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Appendix  A 

TARGET  DETECTION  MODELING  EQUATIONS  (AFTER  BAILEY  AND  MUNDIE) 


VISUAL  DETECTION  EQUATIONS 

H.  H.  Bailey's  visual  target  detection  model  [1]  has  been  somewhat 
abbreviated  and  modified  for  use  with  the  WETTA  algorithms.  The  main 
change  is  that  I have  not  so  far  used  Bailey's  "search  term"  (the  prob- 
ability that  an  obser'er,  searching  the  target  area,  looks  for  a spec- 
ified glimpse  time  in  the  direction  of  the  target  with  his  foveal  vi- 
sion). Therefore,  environmental  effects  on  scene  clutter  (or  congestion) 
have  not  been  dealt  with  explicitly  as  they  would  affect  search  time. 

I have  adapted  Bailey's  "detection  term"  (the  probability  that  a target 
is  seen  given  that  it  is  in  the  observer's  foveal  field  of  view)  and 
"resolution  term"  (the  probability  that  a target  is  adequately  resolved 
to  be  recognized — the  degree  of  recognition  required  can  be  specified) ; 
these  are  ^nd  respectively,  in  my  nomenclature.  The  prob- 

ability of  detection  with  adequate  resolution, 

p = p ♦ p ^ 

DRvis  uvis  Rvis' 


p 

Dvis 

where  "+"  is  for  ^ 1,  otherwise;  is  target-to-back- 

ground  contrast  transferred  through  n atmospheric  layers,  and  the 
threshold  contrast  of  the  observer- target-range  combination.  is  de- 
fined as  the  received  contrast  that  yields  a 50  percent  detection  prob- 
ability, which  Bailey  has  approximated  by  the  hyperbola, 

^°®10  S + 0*5)”^  - 2, 

♦ 

^If  clouds  intervene  between  target  and  sensor,  then  I*DRvis  ~ 
^CFLOS  ' ^Dvis  ‘ ^Rvis’  ^CFLOS  P’^°babllity  of  a cloud-free 

line  of  sight. 


B5BBSS 


where  a is  the  angular  size  (minutes  of  arc)  of  the  small  dimension  of 
the  target  intercepted  at  the  sensor. 

a « 3.A4  M£/R, 

where  M is  magnification » 2.  is  a characteristic  (usually  minimum  pro- 
jected) dimension  of  the  target  (m) , and  R is  range  (km) . The  received 
contrast  transferred  through  n atmospheric  layers 


C «=  C 
n o 


i5i  (fi)i 


where  C is  Inherent  target-to-background  contrast  and  (F.j)^  is  con- 
trast "transmittance"  (transfer  function)  through  atmospheric  layer  1. 


where  (L  )^  is  the  background  reflectance  and  (L  )^  the  target  reflec- 

g .1  t X 

tance  at  the  bottom  of  the  lowest  atmospheric  layer  (i  * 1). 


(^t)i  ■ [i  * (%  • 


where  (SG).,  6,  and  R are  the  estimated  "sky-ground  ratio,  visible 

i i i -1 

extinction  coefficient  (neper  km  ),  and  path  length  (km),  respectively, 
in  atmospheric  layer  1.  For  elevated  layers  (i  > 1),  sky-ground  ratio 
is  calculated  from  an  Inferred  albedo,  a^,  taken  as  being  equivalent  to 
the  background  reflectance  found  at  the  bottom  of  that  layer,  i.e.: 

[(Oi-l  " (^)i-lle 


where 


(^g)i-i " (^t)i-i  * Vi“i-i. 
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The  "resolution  term"  (probability  of  adequate  target  recognition) , 

, -» 2 
I / , 1 \ 

1 - exp ( - 


■Rvis 


a'  0.4F 


- 1 


r 


where  a'  is  required  target  angular  size  (minutes,  at  observer)  for  a 
nearly  1.0  probability  of  seeing  the  target  signal  and  is  a subjec- 
ti''e  "shape  recognition  factor”  ranging  from  1 (spot  recognition;  "some- 
thing is  there")  to  5 (high  confidence  shape  recognition)  by  means  of 
which  "adequate  resolution"  is  specified. 


logfo  ci'  = ilogj^g  (2C^)  f 2) 


-1 


- 0.5. 


If 


a'  0.4F  - 
s 

then  P„  , =0. 

Rvis 

The  above  formulation  is  for  detection  by  the  human  eye.  In  apply- 
ing it  to  TV  seeker  systems,  I have  implicitly  assumed  that  if  the  un- 
aided eye  can  detect  a target,  then  a properly  designed  TV  system  and 
its  operator  can  attain  and  maintain  the  necessary  lock-on  to  the  tar- 
get. 

HR  OETECTIOl)  EQUATIONS 

The  following  set  of  equations  for  HR  target  detection  probabil- 
ity is  taken  from  the  work  of  Bailey  and  Mundie  (In  Lau  et  al.  [2]). 

As  with  the  visual  detection  equations,  the  probability  of  detection 
with  adequate  resolution. 


P = P • P 
DRlr  Dir  Rir’ 


where  is  the  probability  of  seeing  the  target  image  on  the  HR 


As  noted  before,  with  intervening  clouds  P.,„,  = P„„, • P.,, 

® DRir  CFLOS  Dir 


■Rir' 
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display  given  that  the  target  image  is  in  tlie  operator's  foveal  field 

of  view,  and  P is  the  probability  that  the  target  image  is  adequately 
Kir  ^ 

resolved  for  the  required  degree  of  recognition. 


Dir 


exp 


-1 

L 

where  S /N  is  the  resulting  signal- to-noise  ratio  after  image  transfer 
n n 2 

through  n atmospheric  layers.  If 


S /N  <1,  then  P,..  = 0, 

n n — * Dir  ’ 


where  is  target-to-background  temperature  difference  (°C) , is 

atmospheric  8-L2  um  transmittance  through  n atmospheric  layers,  and 

F.,  and  T.„„  are  instrumental  modulation  transfer  function  and  noise 
M AN£ 

equivalent  t(  ■ perature  difference  (°C),  respectively. 


■ l2l  ’!• 


being  the  8-12  urn  transmittance  through  layer  i,  and 


'i  = ^ 


-R. 

i m c a 1 


This  fonri  for  the  noise  dependence  of  Pj)  is  not  the  only  one  used. 

2 

This  equation  for  S/N  accounts  for  detector  noise,  only.  The  for- 
mulation by  Rosell  and  Willson  [6]  accounts  for  other  i.olse  sources, 
and  also  uses  the  Gaussian  form  for  the  equation  ior  Pp.  Comparative 
calculations  (unpublished)  by  T.  F.  Lippiatt  of  The  Rand  Corporation, 
using  bodi  the  above  and  the  Rosell  and  Willson  aquations,  showed  that 
Pp  predicted  by  the  two  methods  could  be  made  nearly  equivalent  by  ad- 
justment of  a single  parameter,  display  gain,  to  a value  near  5. 


^ ' - '-"'0 


where  R is  path  length  through  atmospheric  layer  i and  v , y and  Y 
1 m c a 

are  H^O  molecular  absorption,  H^O  continuum  absorption,  and  aerosol  ex- 
tinction coefficients,  respectively,  in  layer  i. 

The  probability  of  adequate  target  image  resolution, 

Pj^lr  = 1 - expj-  (N^/0.4F^)  - 1 j, 

where  is  the  number  of  HR  display  resolution  lines  across  the  dis- 
played target  and  is  the  shape  recognition  factor  as  previously  de- 
fined; and 

M *4  ..  ■ . 

r 1.19aR  ’ 

where  M,  )!,,  and  R are  the  magnification,  target  dimensions,  and  range 
as  previously  defined,  and  a is  the  angular  size  (mrad)  of  the  display 
resolution  element.  (The  factor  1.19  results  from  applying  the  Kell 
factor,  ^2,  to  the  vertical  resolution  element  and  taking  the  geometric 
mean  of  the  stze  of  the  vertical  and  horizontal  resolution  elements.) 
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